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ABSTRACT 

Theories of turbulent time correlations are applied to compute frequency spectra of sound 
radiated by isotropic turbulence and by turbulent shear flows. The hypothesis that Eule- 
rian time correlations are dominated by the sweeping action of the most energetic scales 
implies that the frequency spectrum of the sound radiated by isotropic turbulence scales 
as u 4 for low frequencies and as u r 4 ^ 3 for high frequencies. The sweeping hypothesis is 
applied to an approximate theory of jet noise. The high frequency noise again scales as 
cj - 4 / 3 . but the low frequency spectrum scales as ui 2 . In comparison, a classical theory of jet 
noise based on dimensional analysis gives uj ~ 2 and u> 2 scaling for these frequency ranges. It 
is shown that the u~ 2 scaling is obtained by simplifying the description of turbulent time 
correlations. An approximate theory of the effect of shear on turbulent time correlations is 
developed and applied to the frequency spectrum of sound radiated by shear turbulence. 
The predicted steepening of the shear dominated spectrum appears to be consistent with 
jet noise measurements. 
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I. Introduction 


The nature of sound radiation by turbulent flow depends on properties of turbulent 
time correlations. For even if the energy contained in motions of any given spatial scale is 
known, say from the Kolmogorov theory, these motions will radiate sound at a frequency 
proportional to their inverse correlation time; this dependence evidently determines the 
frequency distribution of acoustic energy. 

The present work applies theories of turbulent time correlations 1-5 to compute fre- 
quency spectra of sound radiated by isotropic turbulence and by shear turbulence. A 
question which must be answered at the outset is whether Eulerian or Lagrangian time 
correlations are required by this calculation. In the usual picture 6 of quadrupole sound 
sources, sound radiates from every point in an incompressible flow field because infinitesi- 
mal regions surrounding each point change shape without changing volume. In this picture, 
it is the Eulerian time correlations which are relevant. A calculation which treats mov- 
ing fluid particles as acoustic sources must consider the random Doppler shift due to the 
motion. 

Analyses of sound radiation by stationary turbulence have often assumed 7,8 that the 
space and time variables in the space-time correlation function are separable. Thus, for 
homogeneous turbulence, 


< i>j(x,t)vi(x + r, s) > = F(t)R(t/t*) (1) 

where r = t - s denotes time difference and t* denotes an integral time scale. But the 
temporal part of the space-time correlation function is scale dependent and Eq. (1) should 
be replaced by 

< v*(x, t)vi(x + r, s) > = F(r)R(r, r) (2) 

where 

R(r, 0) = 1 (3) 

The goal of this paper is to assess the implications of using the more realistic model Eq. (2) 
in noise calculations. While it may be possible to compute single time quantities including 
the total acoustic power using a simplified formula such as Eq. (1), accurate prediction of 
two-time quantities like the acoustic power spectrum cannot be expected in general. 
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The specific form of Eq. (2) applied here is determined by the sweeping hypothesis for 
Eulerian time correlations, 1,3,4 which is described in more detail in Sect. II. Our previous 
work 9 applied this hypothesis to the calculation of the total acoustic power radiated by 
isotropic turbulence. In this idealized problem, the power spectrum is shown to scale as 
uf~ 4/3 for large frequencies and as u ; 4 for small frequencies. Some numerical simulation 
data consistent with these scalings are described by Lilley. 10 

To demonstrate the importance of time correlations in predicting acoustic radiation 
from inhomogeneous flow, the classical theory of jet noise 11,12 is evaluated using Eqs. (2) 
and (3). This theory applies the Lighthill formula for far-field sound radiation, and as- 
sumes isotropic turbulence statistics. Although this theory is highly simplified, it provides 
a convenient way to compare the consequences of different theories of turbulent time corre- 
lations. The high frequency acoustic power is found to scale as w~ 4 / 3 , like isotropic turbu- 
lence, but the low frequency spectrum scales as cn 2 . In contrast, dimensional arguments 12 
lead instead to u 2 scaling for low frequencies, and uj~ 2 scaling for high frequencies. It is 
shown that these scalings are recovered if the space-time correlation function is given the 
separated form postulated in Eq. (1). The difference in predicted scalings demonstrates 
the role of time correlations in acoustic analysis. 

Although there is some experimental evidence that noise spectra measured at 90 
degrees from the jet axis, for which this type of theory is most accurate, decay somewhat 
more slowly than u;' 2 at high frequencies insofar as they exhibit power law scaling at 
all, these comparison with jet noise data are inconclusive. Quantitative predictions will 
depend on the numerical integration of the linearized Euler equations 13 with the present 
theory used to characterize the sound source. 

Finally, we attempt to evaluate the power spectrum of sound radiated by shear tur- 
bulence. The concern of this analysis is not the so-called “shear noise,” 7 which is a global 
property of the mean flow, but the corrections to the acoustic frequency spectrum due to 
modifications by shear of the space-time correlation function. In contrast, the standard 
theory 7 does not consider these effects, but emphasizes instead the effects of shear on the 
directivity pattern of acoustic radiation. 

Since time correlations in shear turbulence do not seem to have been investigated 
previously, they are computed using Leslie’s perturbation theory 14 for the direct interaction 
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approximation for shear turbulence. 15 As in Leslie’s calculation, the acoustic spectrum 
is expressed as an expansion in powers of the mean velocity gradient. A consequence 
of the two-time properties of turbulence is that the radiated sound can depend on the 
antisymmetric part of the mean velocity gradient alone, unlike the Reynolds stress, for 
which such dependence is ruled out by general invariance arguments. 16 Thus, it may not 
be possible to characterize sound radiation by the Reynolds stress tensor. If the shear 
dominated noise is produced primarily by the near exit region of a jet, its frequency 
spectrum will scale as u )~ 10 / 3 . This prediction compares reasonably well with a recent 
summary 17 of jet noise data. 

II. General properties of acoustic frequency spectra 

According to the Lighthill theory, 6 the frequency spectrum of sound radiated by ho- 
mogeneous turbulence can be written as 


pH — 


(4) 


In Eq. (4), A is Ribner’s 7 geometric factor 

A ijkl = ^ ^ XjXj XfcXjX 

where the vector x t connects the source and observation point, and ^ denotes symmetriza- 
tion by summation or integration over source points symmetric with respect to the statis- 
tical symmetries of the turbulence. The dynamic factor / is 

u ; 4 f 

Iijki(u) = 167r2c5 / dp[Q ik (p,u)*Q j i{-p,u}) + Q i i(p,u})*Q jk (-p,uj)] (5) 

where * denotes convolution with respect to u> and <5tj(p, u;) is the two-time correlation 
function defined by 

< Vi{p,u)Vj(p',Uj') > = < 2 j.,(p,u>)< 5 (c*;-l-n/)< 5 (p + p') 

Angular dependence of p in Eq. (4) arises from angular dependence of both A and Qij. 
Eqs. (4) and (5) are formulas for the acoustic power spectrum per unit mass of fluid; 
accordingly, the mass density p 0 of the fluid does not appear. 
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We recall the arguments which lead to Eqs. (4) and (5). The starting point is 
Lighthill’s inhomogeneous wave equation. The solution of this equation is an integral 
containing space derivatives of time-retarded velocities. The dominant terms in the solu- 
tion are identified as the space derivatives of the time-retarded quantities. These space 
derivatives are equal to time derivatives, which for time stationary conditions contribute 
the factor u; 4 in Eq. (5). 

This argument leads to a fourth order velocity correlation; another approximation 
made in arriving at Eq. (5) is that the familiar quasinormal closure hypothesis is applied 
to close this correlation. An alternative to quasinormality which could potentially alter the 
sound pressure levels is considered in Eef. 18, but the corrections computed are very small. 
Another alternative which could alter the frequency scaling of aerodynamically generated 
sound is anomalous scaling of pressure correlations. 19 At this time, the occurrence of such 
anomalous scaling remains conjectural. 

A simplified theory of sound radiation from shear flows results from a straightforward 
generalization 11 of Eq. (4) to inhomogeneous turbulence, 

p(oj) = J dy Ayfc|/yfci(w,y) (6) 

where y denotes position in flow field and the integration is over the flow region. In this 
case, Eq. (5) is generalized to 

W-.y> = ^ / rf p 

[Qik( P, y- U>) * Qji(- P, y, u) + Quip, y, w) * Q jk (- P, y, w)] (7) 

Since Eqs. (6) and (7) ignore flow-acoustic interaction, 13 their predictive capability is 
limited. These formulas will be used to compare the consequences of different theories 
of turbulent time correlations on noise spectrum predictions, not to obtain quantitative 
predictions. 

III. Frequency spectrum of sound radiated by isotropic turbulence 

For isotropic turbulence, Ribner’s factor is 


A ijkl — dij&kl “A di k 8ji -|- 8u 5k j 


(8) 
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The space-time correlation function of stationary, homogeneous, isotropic turbulence has 
the form 

/ OO 

dr R i0) (p,r)e tWT (9) 

-OO 

where 

Pij( P) = $ij ~ PiPjP~ 2 


The superscript (0) is used in Eq. (9) and subsequently to identify isotropic quantities. 
The single time correlation function 


4nk 2 Q < '°\p) = E(p) ~ 


e 2/3p-5/3 f Qr p > £ o 

E 0 (p) for p < k 0 


( 10 ) 


exhibits Kolmogorov scaling for p > ko; ko is the inverse integral scale which character- 
izes the beginning of the inertial range. The far infrared spectrum Eq is included for 
completeness; however, it will not influence the frequency scaling of the acoustic power 
spectrum. 

The time correlation function has an inertial range similarity form 


R w ( P ,t) = R w (t0( P )) 


( 11 ) 


where the functional form of 6(p) is determined by the temporal decorrelation mechanism. 
For the sweeping hypothesis for Eulerian time correlations, 1,3,4 


0(p) = C D Vp (12) 

where the sweeping velocity V in Eq. (12) is a property of the most energetic scales 
of motion. Thus, under the sweeping hypothesis, temporal decorrelation has a nonlocal 
character since it is determined for motions of any given scale by motions of possibly much 
larger scale. It is therefore possible that whatever contribution the most energetic scales 
may make to the total acoustic power, they always determine the frequency distribution 
of acoustic energy through their effect on the temporal decorrelation of the inertial range 
scales. In free shear flows, it is very likely that the sweeping velocity is a property of 
large-scale coherent structures. In this case, the sweeping velocity will not be predicted by 
turbulence models, which compute either inertial range or mean flow quantities. 
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Substituting Eqs. (8)-(12) in Eq. (5) yields 


Wiv ilpL dT * <0)(e(p)T)V “ 

where the second equality defines R ( ' 0) and Q is the geometric factor 

Q = Aijki j> dS( k) [P ifc (k)P,-,(k) + Pu(k)P jk ( k)] = 64tt 


Since 


/ oo 

dr 7?^(p, r) 2 = R < 

- OO 


the limit of Eq. (13) for small u; gives for the sweeping hypothesis, 


4 /*oo 


£(p) 2 1 


p(u>) = — / d p^- — 
7 r 2 c 5 Jo p 2 Vp 


T his result expresses p as a functional of the exact energy spectrum E(p). A simple 
approximation is obtained by substituting the cutoff Kolmogorov spectrum, Eq. (10) with 


Eo = 0: 


16/3 


For large u;, introduce the change of variables £ = 0(p)/u; in Eq. (13) to obtain 


*0 = 4* f V* TTJ^ km ^ 

tt 2 c & Jo F M;/F) 

£ 4 / 3 ^ 1373 ^- 4/3 


where the Kolmogorov spectrum for p > k 0 from Eq. (10) is again substituted for E(p). 
Summarizing the results of Eqs. (14) and (15), 


P£ ( k~ 16/3 V-h^ 4 for u ~ 0 
c 5 \ e 1 /3p r i3/3 u; -4/3 for u; > Vk Q 


(16) 


In this calculation, the acoustic power radiated by inertial range scales is determined 
through nonlocal sweeping effects by the sweeping velocity F, which is characteristic of 
the most energetic scales. The peak acoustic power occurs at a frequency u> ~ Vk 0 - Some 
preliminary evidence supporting the -4/3 scaling is given by Lilley 10 . 
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It is useful to compare the scalings of Eq. (16) with the scalings which would result 
from the hypothetical relation 

e ~ P a (17) 

with a ^ 1. An example is the straining hypothesis for Eulerian time correlations, 1,4 for 
which a = 2/3. A simple calculation shows that for high frequencies 

p{uj) ~ , j3 = 4 — + a)a~ l (18) 

O 

In particular, for the straining hypothesis, p ~ u~ 7/2 . For any a < 1, the high frequency 
spectrum predicted by Eq. (18) is steeper than the high frequency spectrum predicted by 
Eq. (16). This result can be understood in terms of the properties of time correlations: 
Eq. (12) shows that the correlation time of a given small scale, for which p » 1 in 
suitable units, is shorter under the sweeping hypothesis than under Eq. (17) with a < 1; 
therefore, the sweeping hypothesis predicts that any given small scale radiates sound at 
a higher frequency. This places more acoustic energy in higher frequencies and therefore 
makes the spectrum more shallow. 


IV. Approximate theory of jet noise 

In this section, the approximate scaling theory of jet noise 11,12 will be modified by 
introducing the time correlation function defined by Eqs. (11) and (12). In this theory, the 
sound radiated by an inhomogeneous flow is given by Eq. (7) in which all quantities are 
functions of position. The calculation will assume that the turbulence is locally isotropic. 
We also use the result 11 that the mean motion causes a Doppler shift of the acoustic 
spectrum as a whole, but does not alter the frequency scaling of the spectrum. 

The theory uses similarity solutions to describe the flow field and suppresses all vari- 
ation in cross-stream planes by averaging. Thus, all flow variables f(x,y,z) are replaced 
by average values 

I dxdz f(x,y,z) 

where A(y) denotes the jet area at downstream distance y. The theory further assumes 
that there are two similarity regimes for a jet: the exit and the fully developed region. To 
simplify the notation, these regions will be taken to extend over 0 < y < L and L < y < oc 
respectively where y is the axial coordinate; the transition region between these similarity 
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regimes will be ignored. L can be assumed to equal about 6.0£> where D denotes the exit 
diameter of the jet. 

A. The fully developed region 

The classical similarity assumption for a fully developed jet is that the total momen- 
tum flux 

M =»S dxdzU(x ’ y ' z? 

is constant downstream and that all cross-sectional averages scale with A4 and y only. In 
terms of the exit velocity U and diameter D. 

M ~ U 2 D 2 (19) 

In particular, this hypothesis gives the well-known scaling for average velocity 

U(y)^M 1 / 2 y - 1 

and also determines the downstream variation of cross-section averages of turbulence quan- 
tities; thus, 

i(y) ~ A4 3/ V 4 (20) 

It follows from Eq. (20) that the cross-section averaged spectrum satisfies 

Q (0 \k,y)~ My-* /z k~ ll/z (21) 

In Eqs. (20) and (21), ~ denotes equality up to a universal constant; all dimensional 
parameters appear on the right hand side. In principle, these constants can be found from 
measurements. 

Denote by ko(y) the integral scale at downstream location y. Then evidently. ko(y) = 
ajy where a is a universal constant proportional to the inverse spread rate of the jet. The 
total energy is found by integrating Eq. (21) over all k > a/y and leads to 

K(y)~My ~ 2 

Introduce following Lilley 11 

p(y, w) = J dxdz p(u », x, y, z) = A(y)p(y , w) 
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Under the sweeping hypothesis, the cross-section average of the frequency convolution in 
Eq. (7) has the similarity form 

A(y) J dxdz R{0) ( U / V P) * R {0) (“/Vp) = ~ ^ 

The subsitution of either an eddy turnover time like Kje or the inverse mean strain rate 
for the sweeping velocity will alter the proportionality constant in Eq. (22). It is likely that 
use of either of these time scales in calculations will require adjustment of this constant 
for each application. 

Now substitute the results of Eqs. (21) and (22) in Eq. (7), then integrate over 
wavenumber p from a jy to oo. The result of this calculation is 

p(y, U ) ~ c-VM'Vy - 7 ' 3 jf” dp 

= (23) 


In Eq. (23), G is defined by 


G( 


<■ *>r 


H 


2 A^f-1/2^-1 


ury 4 M 


<K C a/3 F(0 


M}! 2 

Note that by integrating over u, we obtain 

/ OO 

du p(y.u) ~ c- 5 jW 25 /V 23 / 3 (!/ 2 .M- 1 / 2 r 1/3 ~ c-*M \- 7 

-oo 

which is the y ~ 7 law for power radiated by each cross-section. 12 
For small arguments, 


G(u>y 2 M~ l/2 ) ~ (u>y 2 M ~ 1/2 ) 16/3 


(24) 


For large arguments, if F decays exponentially or faster at infinity, then 

G(uiy 2 M~ 1 ^ 2 ) ~ constant (25) 

The power spectrum p(u) is evaluated by integrating Eq. (23) with respect to y ; since the 
integral converges at y = oo, it can be extended to oo for analytical convenience, and the 

result is oo 2 

p(u) ~ c- 5 u,-'W 5 ' 6 r dy ! r a/ 3 G(j^) (26) 
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In view of Eq. (25), this integral converges in the limit u > — > oo; therefore, in this limit 


p(u) ~ w -4/3_ M 25/6 L - 20 /3 


(27) 


To evaluate the limit of Eq. (26) when u> is small, it is convenient to introduce the 
change of variable Q 2 = ujy 2 Ai~ 1 ^ 2 . The result is 

rOO 

PM - c- 5 u 2 M 5/2 / dC G(C 2 )C ~ 23/3 

Ju'/iLM - 1 / 4 

In view of Eq. (24), G((j 2 ) ~ C 32 ^ 3 for small (, therefore, this integral converges for small 
arguments and when u — > 0, 

p(uj) ~ a; 2 A4 5 ^ 2 (28) 


Summarizing Eqs. (27) and (28) 


p(uj) ~ c 


t/ 5 T> 5 u 2 u ~ 0 

^ 25 / 3 ^ 573 ^- 4/3 u ^ 


(29) 


where Eq. (19) has been substituted in Eqs. (27) and (28) and the result L ~ D is used 
to eliminate L. 


B. The exit region 

The exit region of a jet has the scaling properties of an axisymmetric mixing layer: 
there is a constant velocity scale, and the area of the mixing region increases linearly 
downstream. Thus, 



U(y ) ~ y° 

(30) 


A(y) = J dxdz ~ yD 

(31) 

Then the scalings of Eqs. 

(20) and (21) are replaced by 



U z y~ l 

(32) 


Q m (p,y) ~ u 2 y- 2 / 3 p - ll/3 

(33) 

and 

r dxdz R m (u>/Vp) • R m (w/Vp) ~ 

Up Up 


mi 

(34) 
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Substituting Eqs. (30)-(34) in Eq. (7), 


p{y, u) ~ u ;- 4/3 c- 5 y- 1/3 f/ 25/3 M(o; 2 /t/- 1 ) 


(35) 


where 


rujy/U 

H(uy/U) = / dli £ 13/3 F( 0 

Jo 

Integrating Eq. (35) over u ; leads to 

J du p(y , w) ~ y~ 1 ^ 3 U 25 ^ 3 D(yU~ 1 ) 1 ^ 3 ~ U S D 

which is consistent with the y° law for power contributed by each downstream section. 12 
Integrating Eq. (35) with respect to y, 

r u>L/U 


(36) 


p(uj) 


rUJLif u 

Du~ 2 U 9 c ~ 5 / C _1/3 ^(0 

JO 


(37) 


In view of Eq. (36), H(Q is constant in the limit £ -4 oo. The integral in Eq. (37) 
therefore diverges as £ 2//3 for large arguments; accordingly, for w -> oo 


p(w) ~ U!- 2 C- 5 U 9 D{uL/U) 2/3 ~ ^-4/3,725/3^5/3 


(38) 


Note that the divergence influences the scaling exponent, which therefore cannot be ob- 
tained by dimensional analysis alone. For small arguments, H(() ~ £ 16/3 , therefore at low 
frequencies, 

pH ~ Du>~ 2 U 9 (ujL/U) 6 = U) 4 U 3 DL 6 (39) 


Summarizing the results of Eqs. (38) and (39), 

/ \ -5 / t/ 3 £> 7 w 4 w ~ 0 

1^25/37)5/3^-4/3 w ^ 


(40) 


Comparing Eqs. (40) and (29), it is evident that the oo 2 contributed by the fully developed 
region is dominant at low frequencies. Combining the contributions from the exit and fully 
developed regions, 

f x -5 / U 5 D 5 uj 2 u> ~ 0 

P(w)~C | ^25/37,5/3^-4/3 w ^ ( ) 

The conclusion that the fully developed region dominates the low frequency regime will 
be seen to agree with the classical theory of jet noise. But unlike the classical theory, the 
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high frequency region in Eq. (41) contains inertial range contributions from both the fully 
developed and the exit regions of the jet. 

V. Comparison with the classical theory of jet noise 

An elementary scaling theory of jet noise 12 results from the additional hypothesis that 
each cross-section of the jet emits sound at one dominant frequency. Dimensional analysis 
based on the scaling laws Eqs. (19) and (30)-(31) at once implies that the contribution of 
the exit region to the acoustic frequency spectrum scales as u >~ 2 and that the contribution 
of the fully developed region scales as u/ 2 ; analysis of Lighthill’s integrals as in Sect. IV is 
not required. These are the well-known +2 and -2 scaling laws. 

A refinement of this analysis is due to Lilley, 11 who argues that although the contri- 
bution to the frequency spectrum from each cross-section extends over a finite range of 
frequencies, the observed spectrum is the envelope of these contributions, and again scales 
as u~ 2 . 

These results will be rederived by suitably simplifying the time correlation functions. 
This calculation based on Light hill's integrals will also provide a slight refinement since it 
provides formulas valid for all frequencies; it will be possible to show explicitly that the 
contribution of the fully developed region is vanishingly small for large uj and that the 
contribution of the exit region scales as u/ 4 for small u>. 

The fundamental idea of the classical theory can be incorporated in the present anal- 
ysis by suppressing the wavenumber dependence of the time correlation functions; thus, 
Eq. (22) for the fully developed region is replaced by 

^ J dxiz R«V/V» . R^MVp) = (41) 

and Eq. (34) for the exit region is replaced by 

-1- J dxdz RV>Hu,/Vp) , R™(<j/Vp) = J}F(^) (42) 

Thus, we effectively assume the separation of variables form Eq. (1) for the space time 
correlation function. 

Repeating the calculation leading to Eq. (23), 

rOC 

p{u,y) ~ u 4 c~ 5 M 3/2 y~ 4/3 / p~ 16/3 F(ujy 2 M~ 1/2 ) 

•'i/y 

= u! 4 c~ 5 M 3/2 y 3 F(ujy 2 M~ 1 ^ 2 ) (43) 
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Integrating Eq. (43) over y leads to 


rOO 

p{u>) ~ to 4 c~ 5 Ai 3 ^ 2 J dy y 3 F(ujy 2 Ai~ 1 ^ 2 ) 

/ OO 

d(, £F(£) 

The integral in Eq. (44) converges when u = 0; accordingly, for small uj, 


( 44 ) 


p(u>) ~ L 0* 


(45) 


in agreement with the ”+2” law. For large uj, assuming rapid decay of the time correlation 
function F for large arguments, instead 


p(u) ~ 0 


(46) 


In the near exit region, Eq. (42) leads to 

p{tu, y) ~ uj 4 c~ 5 U 3 y b DF(u>y/U) (47) 

Integrating Eq. (47) over y. 

p(u) = u> 4 c~ b U 3 D [ dy y 5 F(ujy/U) 

Jo 

rUjL/U 

= uj 4 c~ 5 U 3 D(U/uj ) 6 / d£ Z 5 F{0 (48) 

Jo 

Again assuming sufficiently rapid decay of the time correlation function F at infinity, the 
integral converges asw->oo leading to 

p(w) ~ U 9 Dlo~ 2 (49) 

for large uj, in agreement with the —2 law. For small c a, the integral is of the order ( ujL/U ) 6 ; 
consequently, 

p(uj) ~ uj 4 (50) 

for small uj. Combining Eqs. (45), (46), (49), and (50), we recover the classical scalings 

U '~° (51) 

I UJ UJ ~ OO 
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VI. Frequency spectrum of sound radiated by shear turbulence 

This section investigates the effects of shear on turbulent time correlations and the 
consequences of these effects for the acoustic power spectrum. The emphasis on time 
correlation effects should be stressed; previous investigations of sound radiation from shear 
turbulence have analyzed instead the effects of shear on directivity patterns, whether these 
effects are due to local anisotropy, 21 or due to mean flow effects. 7 

The effects of shear on time correlations will be evaluated by applying Leslie’s per- 
turbative theory of shear turbulence. 14 In this theory, the shear is treated as a formally 
small perturbation of a background state of isotropic turbulence; one can compare also the 
two-scale theory of Yoshizawa. 22 The result of the analysis is a perturbation series of the 
form 

<Mk.l. s) = <4 0) (k, t,s) + ») + <?i 2| (k, «,*)••• (52) 

where is of order p in the mean velocity gradient. The single time correlation was 
computed to first order by Leslie; 14 for the corresponding calculation to second order, see 
Ref. 23. Corresponding to Eq. (52) is the series for the acoustic spectrum 

p(ui) = + p^(w) + P^W) H 

where p (0) ( u) is the spectrum of isotropic turbulence. 

Substituting the perturbation series Eq. (52) in the expression for acoustic radiation, 

4 p rOC 

dp L dTe ^ x 

{oji’fp, p, t) + Qj'lp. t)Q { °J(- p, t) 

+ Q'“’(p. r)<$>(-p, r) + Q<, 0) (p, t)Q£(-p, t)} (53) 


The index symmetry of Ribner’s factor implies that Eq. (53) can be rewritten as 

4 r roo -i 

P (1) M= 7^2~s Xjkl dp drei “ T x Q\k(P’ T )Qfl } (-Pi T ) 

1D7r J J ~ co (■ ijkl ) 

where denotes summation over all index permutations. The factor of 1/6 arises 

from the four terms in Eq. (53) and the 24 permutations in YL(ijki)- ^ substituting the 
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expression for r) from the Appendix, a simple calculation shows that 

f dS{p) Q\ 1 k l (p,T)Qf l \p,T) 

(ijki) J 

90 

= -Jg [R + {p, t) + R~{p, -t)]R { 0 ) {t)Q(p ) 2 ^2 s i} 6 ki (54) 

(ijkl) 

The time correlation functions R + and R~ in Eq. (54) are defined in the Appendix. 

Assuming forms for the time correlation functions R( 0 \R + ,R~ consistent with the 
sweeping hypothesis leads to 

pl»(w) (55) 

where 

^ (Sij) = Aijkl ^2 Sij^kl 

(ijkl) 

is a linear invariant of the rate of strain tensor, invariance being understood with respect 
to the statistical symmetries of the flow. 

The perturbation series can be continued in principle to higher order terms in the mean 
velocity gradient. Because the results are very lengthy, let us attempt to generalize Eq. (55) 
heuristically. It is noted in the Appendix that whereas the single time correlation function 
<2^*(k, t, t) can only depend on the symmetric part Sij of the mean velocity gradient, the 
two time quantity Q\j\k, t, s ) depends in general both on Sij and the antisymmetric part 

w _i,m _su i) 

"~2 y dx, dXi’ 

Although the index symmetry of Ribner’s factor prevents dependence of p( x )(u>) 

on Wij, p^{w) certainly can depend on both isotropic invariants SijSij and W,jWij. 
Dependence of the Reynolds stress itself on WijWij can be ruled out on general invari- 
ance grounds. 16 Thus, it is possible that acoustic radiation by shear turbulence cannot be 
parametrized in terms of the Reynolds stress. From the viewpoint of this paper, such de- 
pendence would be problematic in any case since the Reynolds stresses are not associated 
with any temporal decorrelation mechanism. We conclude that 

pW{u>)~i 2 {S ij ,W ij )u- l0 f 3 (56) 
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The approximate calculation of jet noise of Sect. IV can be repeated using the second 
order perturbation series Eq. (52) to describe the sound source. Like the result for isotropic 
turbulence statistics, the results to first and second order will coincide with Eqs. (55) and 
(56). A turbulent jet has the symmetry type of axisymmetric shear flow in which Sij is of 
the form 


Sn Si 2 o 
Sl2 S22 S23 

0 S23 Sn 


(57) 


so that 2Sn + S22 = 0. The largest velocity gradient in a turbulent jet is the shear 
dll 2 /dx\ in the near exit region. If this term dominates noise production, then in view of 
the structure of Eq. (57), I 1 (Sij) ~ S 22 ~ 0 but I 2 (S lj ) ~ S? 2 and the shear contribution 
to the acoustic power spectrum will be given by Eq. (56) with u~ 10 ^ 3 scaling. 


These results for the spectra of sound radiated by isotropic turbulence and by shear 
turbulence combine to determine the complete pattern of acoustic radiation through the 
known directivities 7 of sound radiated by isotropic and axisymmetric turbulence. At 90 
degrees to the jet axis, the shear contribution is entirely suppressed, and only the isotropic 
contribution, which scales theoretically as w -4 ^ 3 , should be observed. The effect of shear 
is greatest for nearly on-axis measurements; acoustic spectra measured in these directions 
should scale as ui ~ 10 / 3 . This very simple description of jet noise spectra is qualitatively 
consistent with the comprehensive survey of jet noise data of Tam et al, 17 who find that 
nearly on-axis spectra scale as u~ a with a ~ 2.8. Unlike the 90 degree spectra, the nearly 
on-axis spectra exhibit unambiguous power law scaling. Of course, this comparison is 
only qualitative; more precise comparisons will depend on numerical simulations now in 
progress in which the present theory is used to characterize the sound sources. 


We stress that in this account of jet noise, the “two components of turbulent mixing 
noise” described by Tam et al are consequences of the isotropic and strain dependent parts 
of the space-time correlation function of shear turbulence. To this extent, our theory 
remains a quadrupole theory, although coherent structures, which are prominent in the 
theory developed by Tam et al, enter crucially as the mechanism of temporal decorrelation. 
We propose that the failures of quadrupole theories described by Tam et al might be 
attributed not to inadequacies of these theories as such, but to the incorrect treatment of 
turbulent time correlations with which they have been applied. 
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VII. Conclusions 


Proper parametrization of turbulent time correlations is crucial to calculation of power 
spectra of acoustic radiation from turbulent flows. The present suggestion, based on the 
sweeping hypothesis for Eulerian time correlations has been applied to the simplified clas- 
sical theory of jet noise. 

This theory could be incorporated in computer simulations in two ways. First, it 
permits stochastic synthesis of noise sources for linearized Euler calculations. 24 A sec- 
ond possibility is to incorporate these models in large eddy simulations as noise sources 
corresponding to the subgrid scales. 

Finally, we note that our results favor the computational characterization of noise 
sources by methods capable of capturing the properties of coherent flow structures like 
large eddy simulation because the sweeping velocity depends on coherent flow structures. 


Appendix. Leslie’s theory of shear turbulence 


We summarize the application of Leslie’s perturbation theory of shear turbulence to 
the two time correlation function in shear turbulence. 

We will consider a time stationary shear turbulence. This requires that energy be 
removed from large scales by some external agency; otherwise, as Kraichnan predicted 14 
and experiments 25 later confirmed, the shear will continually create new scales; in the log 
layer, the external agency is turbulent diffusion. 


Write the equation for velocity fluctuations in shear turbulence as 

Uj(k, x EtTTifi (k) / dpdq u TO (p, t)u n (q, t) -(- «Sirn(k, i)u77j(k, t) 

^ */ k=p-fq 

where 


d_ 

dt 


Pimn — k m Pi n ( k) T /j n F > j n (k) 

Pij{k) = Sij - kikjk ~ 2 
and the mean shear operator is 

Sim(t) = -f —(t) + 2k- 2 k i k p ^(t) + 8 im k s ^{t)-^- 

OX m (/ X771 (yOOf* (JrC r 

Leslie writes the solution of Eq. (58) as a perturbation series 


(58) 


Ui(k,<) = ttj-°^(k,t) + u (1) (k, t) -I 
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where u {n) is of order n in the strain rate and the zero order field is isotropic turbulence. 
Leslie’s analysis assumes that the response function 1 ’ 14 ’ 15 G (0) of isotropic turbulence is 
known. If the shear is treated as a perturbation, then properties of the response function 
imply that the first order velocity field is 

u\ l) (Kt) = [ ds G {0) (k,t,s ) (k,s) u£> (k, s) (59) 

Jo 

The two time correlation function in weakly sheared turbulence can be evaluated from 
Eq. (59). To first order, 

Q[j ) (k,t,s)<5(k + k') = < u\ l) {k,t)uf ] {k',s) + u- 0) (k, ^^(k', s) > 

Then 

Q<}>(k, t. s) = pV( k, t,s) 4- ng>(k,t, s) + T^(k, t, s) (60) 


where 


dUi 

dx r 


P^\k, t, s) = - j dr G(k,t,r)^^-Q m j(k,r.s) - j dr G{k, s,r)£^-Q mi (k,r. t) 


f 

Jo 


dUj 
' dx m 


n[j } (k, t,s) — 2 [ dr G{k,t,r)hik p k' 

J Jo 


' 2 f^Qmj(k, r, s ) 

dxjn 


f s dU 

4-2 / dr G(k, s,r)kjk p k~ 2 —-£-Q m i(k,r 1 t) 

Jo UXm 


Tp\k,t,s) = j dr G(k, t, r ) < uf ] {-k, s)k p ^ J^u\ 0) (k,r) > 


p dx r dk r 

dU n d 


+ j dr G(k. s, r) < uf\-k,t)k p -^-^uf\k,r) > 

J 0 - - 7 ■' 7 


The time stationary forms for P-p and II are 


*? 


P^(k.r) = ~R-{k. T )^p-Q(k)P mj (k) - R + (k,r)^J-Q(k)P mi (k) 


dl h 


dUj 


v "’’ / " v 7 dx m 

n^(k,r) = 2 R-(PT)^k t k p k- 2 Q(k)P mj (k ) 


dx r 


dU P 

’ dx m 
dU„ 


+ 2R + (k.T)—lk j k p k- 2 Q(k)P rni (k) 
ox m 


(61) 


(62) 


where 



dr G(k,r)R (0 ^{k,r - r) 
dr G(k , r)f?^°^(fc, r + r) 
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where is the time correlation function of isotropic turbulence. 
To evaluate correlations of the form define 

Aijr (k, t , s) — U{ (k, t) — ^ 

so that 


rt r\TT 

T^ 1 (k, t, s) = - I dr G(k,t,r)kp-~A jir (~k,s, ‘ l 

f s dU 

-J dr G{k,s,r)k p -^A ijr (-k,t,r) 


Note the properties 


(63) 


(k) — 0 

kjAijri k) = -Qtr(k) = - P ir {k)Q(k ) 

Aijr{ k) — Aij r (k) 

r\ rj 

A ijr (k) + A jir ( k) = Qj-Qiji k) = Qj-PijOf)Q{k) (64) 

and the general form, required by isotropy, 

Aij r (k) = k' 1 A(k)eij r + k~ 2 B(k)kiSj r + k~ 2 C(k)kj5i r 

+ k~ 2 D(k)k r 5ij + k~ 4 E(k)k t k 3 k r (65) 


It follows from substituting the general form Eq. (65) into the conditions Eq. (64) that 
the functions A,- - ■ E are uniquely determined so that 

< Ui(k,t)-Uj(-k,s) >= 

Cs rCj-* 

- Q(k,t, s)k~ 2 kjP ir ( k) + t, s)k~ l k r Pij(k) (66) 

Substituting Eq. (66) in Eq. (63) and using Eqs. (60)-(62), we obtain the time stationary 
result 


QSj } (k.r) = Q(k){-R-(k,T)^-P mj ( k) - R+(k,r)^-P mi ( k) 

OXfn OXjfi 

+ 2R~(k, T)^kik p k- 2 P mj (k) + 2R+(k, T)^kjk p k- 2 P mi {k) 

OX m GX m 

- R-{k : T)k- 2 k p kiP jr { k)^ - R + {k 7 T)k- 2 k p k 3 P ir {k)^-} 

+ ±[R-{k,r) + R+(k 7 T)}k- l k r k p P t3 { k)|^ 
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where 


r°° d 

R~(k,r) = I dr G{k,r)—Q{k,r - t) 

r°° d 

R + (k, t) = J dr G(k,r)—Q(k,r + t) 

The symmetry properties of the result in Eq. (67) are interesting. Note the result of 
spherical integration 

I dS(k) Q*‘ , (k,r) = - TQ(t){fi+(i,T) + fi-(t,T)}Sij - i<3(*){B + (t,r) + fi-(*,r)}Sy 

- ±Q(k){R + (k,T) - R-(k,r)}W„ ( 68 ) 

6 

and the single time reduction 

I dS(k) Qlj’(k) = -^Q(k)fi + (fc,0)5y - ±Q(k)R+(k,0)S„ 

Thus, whereas the single time integrated correlation function depends only on Sij, a con- 
sequence of index symmetries alone, the corresponding two time quantity depends on both 
Sij and Wij. This dependence is possible because the odd parity of the function of time 
difference which multiplies Wy in Eq. (68) insures the symmetry relation 

Qij (k, t ) = <5jt(k, — t) 


satisfied by the two-point two-time correlation function. 
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